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Abstract 



Third-order bound-charge electronic nonhnearities of Si nanocrystals (NCs) embed- 
ded in a wide band-gap matrix representing silica are theoretically studied using an 

C^ ' atomistic pseudopotential approach. Nonlinear refractive index, two-photon absorp- 

tion and optical switching parameter are examined from small clusters to NCs up 
i-A , to a size of 3 nm. Compared to bulk values, Si NCs show higher third-order optical 

^ \ nonlinearities and much wider two-photon absorption threshold which gives rise to 

^ ■ enhancement in the optical switching parameter. 
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^ . Both the subjects of Si nanocrystals (NCs) and nonlinear optics are currently 

O ■ very active because of their well established applications, such as those in 

light emitting diodes, lasers, solar cells, interferometers, optical switches, op- 
ts^ ; tical data storage elements, and other photonic devices [1]. One group of very 
\^ • important optical nonlinearities is the third-order nonlinearities which involve 



nonlinear refraction coefficient or optical Kerr index n2 and two-photon ab- 
sorption coefficient (3. These nonlinearities are crucial in all-optical switching 
and sensor protection applications |2j as well as in the up-conversion of the 
sub-band-gap light for the possible solar cell applications [3]. In these respects, 
a clear understanding of the third-order nonlinearities in Si NCs would play 
an important role in such applications. 

Si NCs with controllable sizes have electronic structure largely affected by 
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the quantum and dielectric confinements. Hence, they are expected to have 
markedly different nonlinear optical properties with respect to bulk Si which 
itself already displays significant third-order optical nonlinearities [4,5]. In 
fact, recent experimental reports show that Si NCs have promising nonlin- 
ear optical properties and device applications [GfTfS] . Therefore, an in-depth 
knowledge of Si NCs' nonlinear parameters is essential for various nonlinear 
optics applications. Unfortunately, there is neither a comprehensive theoretical 
work nor an experimental study on the full wavelength and size dependence 
of the third-order optical nonlinearities of Si NCs. Among the few available 
experimental studies, we should mention the work of Prakash et al. in which 
n2 and (3 were measured [6]. However these measurements were performed at 
a single wavelength. For a more comprehensive understanding, a rigorous the- 
oretical work is highly required that can also unambiguously extract the NC 
size effects which is not precise in the experimental studies due to limitations 
in size control in embedded NCs. 

In this paper, our aim is to present such a theoretical assessment of the third- 
order nonlinearities in Si NCs resolving the size scaling and the full wavelength 
dependence up to UV region. In this respect, we indiscriminately consider 
both nonresonant and resonant nonlinearities. Furthermore, unlike most of 
the previous theoretical studies on linear optical properties which considered 
hydrogen-passivated NCs, we deal with Si NCs embedded in a wide band-gap 
matrices representing silica which is the most common and functional choice 
in the actual structures pj . The source of optical nonlinearity in this work is 
the bound (confined) electronic charge of the NC valence electrons filling up 
to the highest occupied molecular orbital. Especially for ultrafast applications 
in the transparency region this is the dominant contribution [9]. 

The characterization of the third-order nonlinear susceptibilities up to a pho- 
ton energy of 4 eV brings a challenge for the electronic structure. For this pur- 
pose we resort to the so-called linear combination of bulk bands basis within 
the pseudopotential framework [10]. This can handle thousands-of-atom sys- 
tems both with sufficient accuracy and efficiency over a large energy window 
which becomes a major asset for the third-order nonlinear optical susceptibil- 
ities with very demanding computational costs of their own. Further credence 
for this particular computational framework is recently established from suc- 
cessful applications of quantum phenomena taking place over several eV en- 
ergy range, such as the excited-state absorption [TT] and the Auger and carrier 
multiplication [12] in embedded Si and Ge NCs. As our primary interest is on 
the nonlinear optical properties, we refer to our previous work for further de- 
tails on the electronic structure [TT]. However, some information regarding the 
embedding host matrix would be in order. In real applications, the common 
choice is silica [1]. On the other hand, to avoid complications arising from its 
chemical and structural compositions, we prefer to replace it with an artificial 
medium having the same crystal structure as silicon but possessing the band 



alignment and refractive index compatible with Si02. Both of these are crucial 
for the accurate representation of the quantum and dielectric confinement of 
the actual system. 

In this work, the electromagnetic interaction term in the Hamiltonian is taken 
as — er ■ E, in other words, the length gauge is used. The third-order op- 
tical nonlinearity expressions based on the length gauge have proved to be 
successful in atomic-like systems [2]. Therefore, we use the expressions based 
on the length gauge whose detailed forms can be found in Ref. [13]. For em- 
bedded NCs an important parameter is the volume filling factor fy of the 
NCs /j, = Vjyc/ysc where Vjyc and Vsc are the volumes of the NC and the 
(computational) supercell, respectively. Throughout this work, for the sake of 
generality the nonlinear properties of the embedded NC systems are calculated 
at the unity filling factor which is indicated by an overbar 
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where the subscripts {a, b, c, d} refer to cartesian indices, u^, ojp, and oja are 
the input frequencies, oj'i = oj^ + up + Ua. Our results can trivially be con- 
verted to any specific filling factor. As we consider only spherical NCs, we 
set the cartesian tensor indices as {a, 6, c, d}={l, 1, 1, 1} and suppress these 
subscripts for convenience. Furthermore, we use a Lorentzian energy broad- 
ening parameter of 100 meV at full width. The value refers to the typical 
room temperature photoluminescence (hence, interband) linewidth of a single 
Si NC embedded in Si02 as considered in this work [Hj. This parameter not 
only broadens the resonances but also introduces the band tailing effects [15] 
which becomes especially important in the transparency regions. 

At high illumination intensities, third-order changes in the refractive index and 
the absorption are observed due to the virtual and real excitations of the bound 
charges. Accounting for these effects, the refractive index and the absorption 
become, respectively, n = no + n2l and a = oq + 131, where uq is the linear re- 
fractive index, a^ is the linear absorption coefficient, and / is the light intensity. 
n2 and (3 are proportional to Rex*-^^(— ^; tu, —OJ., uj) and \tolx}~^\—uj; cu, —uj, cu), 
respectively [16]. Furthermore, the degenerate two-photon absorption cross 
section [2] at unity filling factor is given by 'a^'^\uj) = a^^\uj) / f^; a^^\uj) and 
P are related to each other through /5 = 2hLU'a^^\Lu) . 

In the case of NCs, one should take into account the local field effects (LFEs) as 
in composite materials the dielectric mismatch between the constituents may 
lead to remarkably different optical properties [17] . Incorporation of the LFEs 
into calculations is not a trivial task. For structures of the so-called Maxwell- 
Garnett geometry ^7\, the LFEs yield the following correction factor for the 



third-order nonlinear optical properties [T7] 



3e^ 



^Nc + 2e/i 



?>ei 



(■Nc + 2e/i 



where e/i and e^c are the dielectric functions of the host matrix and the NC, 
respectively. Note that in this equation it is assumed that the inclusions are 
non-interacting spheres and the host matrix does not show any significant 
nonlinearity [18]. In our implementation we use a static local field correction, 
otherwise the correction factor spuriously causes negative absorption regions 
at high energies. Further theoretical work is much needed in this direction. 
Regarding the other shortcomings of our model, we would like to mention 
that our treatment does not include the excitonic, strain, thermal and free- 
carrier effects, which may form the basis of possible extensions of this work. 
Nevertheless, we can argue that in our context the lack of these effects will not 
have qualitative consequences. First of all regarding the strain, a very recent 
and realistic atomistic modeling of this system has revealed that the inner core 
of the NCs where most wave function localization occurs remains unstrained 
[19j. As for the thermal effects, in this work they are indirectly accounted 
through the Lorentzian broadening parameter which is predominantly caused 
by the low energy acoustic phonons [20]. On the other hand, the free-carrier, 
phonon-assisted or other thermal nonlinearities are totally left out of the scope 
of our treatment which focuses on the bound-charge nonlinearities taking place 
in much faster time scale. Finally, the excitonic effects may introduce new 
features to the spectra, however, at room temperature we do not expect them 
to markedly stand out. 

We consider four different diameters, D =1.41, 1.64, 2.16 and 3 nm. Their 
energy gap Eq as determined by the separation between the lowest unoccu- 
pied molecular orbital (LUMO) and the highest occupied molecular orbital 
(HOMO) energies show the expected quantum size effect [11]. The n2 is plot- 
ted in Fig. [1] which increases with the decreasing NC size for all frequencies. 
The smallest diameter gives us the largest n2. When compared to the n2 of 
bulk Si in this energy interval (~ 10"^"^ cm^/W) [lf5|2T] . our calculated n2 
is enhanced as much as (~ 10^/„) for the largest NC. For Si NCs having a 
diameter of a few nanometers, Prakash et al. [6] have obtained ^2 of the order 
of ~ 10~^^ cm^/W which in order of magnitude agrees with our results when 
a typical f^ is assumed for their samples. 

We have plotted (3 against the photon energy in Fig. [21 Peaks at high en- 
ergies are dominant in the spectrum and (5 decreases with the growing NC 
volume. The obtained (3 is about 10^/^ cm/GW for the largest NC at around 
1 eV. When compared to the experimental bulk value (1.5 — 2.0 cm/GW mea- 
sured at around 1 eV) [3], our calculated (3 is enhanced about 10^/^, times. 
Prakash et al [6J have observed /3 to be between (10^ — 10^ cm/GW) at 1.53 eV 
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Fig. 1. n2 as a function of the photon energy for different NC sizes. 
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Fig. 2. /3 as a function of the photon energy for different NC sizes. 
which is close to our values provided that /^ is taken into account. Moreover, 



\x^^\—(^; (^, —uj, uj)\ at 1.53 eV in our results (not shown) is enhanced with the 
decreasing volume. This is in full agreement with the findings of Prakash et al 
[6j. However, the result of the composite material was given in this experiment 
rather than the contributions of each constituent, that is, the host matrix and 
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Fig. 3. Optical switching parameter n2//?A as a function of the photon energy for 
different NC sizes. The symbols refer to experimental values for bulk Si indicated 
by diamond [21] , square [21] , circle |4j and triangle ^ and NC Si indicated by cross 

m- 

Si NCs. 



We should note that (3 is nonzero down to static values due to band tailing 
as mentioned above. Another interesting observation is that the two-photon 
absorption threshold is distinctly beyond the half band gap value, which be- 
comes more prominent as the NC size increases. This can be explained mainly 
as the legacy of the NC core medium, silicon which is an indirect band-gap 
semiconductor. Hence, the HOMO-LUMO dipole transition is very weak es- 
pecially for relatively large NCs. We think that this is the essence of what is 
observed also for the two-photon absorption. As the NC size gets smaller, the 
HOMO-LUMO energy gap approaches to the direct band-gap of bulk silicon, 
while the HOMO-LUMO dipole transition becomes more effective. 

In optical switching systems ?t,2 has an important role f2]. A good optical 
switching device should possess the condition An > Cg^aA where An is the 
change in the refractive index and Csw is a constant of the order of unity, the 
exact value of which depends on the switching system [H] . For the energies be- 
low the band gap, one has Aa = j3I and An = 77-2/. In this case, the condition 
reduces to n2//3A > Csw We have plotted this optical swicthing parameter 
n2/ 13\ as a function of the photon energy in Fig. [31 Their general behav- 
ior resembles the results of Khurgin and Li for the third-order intersubband 
nonliearities in quantum wells 



The condition gives us values already exceeding unity and it brings about an 
immediate peak at around 1 eV for each diameter. We should note that this 
is not the case for bulk Si which shows a monotonic behavior, as computed 
by Dinu for the phonon-assisted two-photon absorption [23j- It can be also 
observed from this figure that the peak positions are red-shifted with the 
increasing NC volume. Notably, the switching condition gives higher values 
as the NC gets smaller, but it converges to a specific value and does not get 
enhanced further. In the same figure some experimental data for the bulk 
Si piH^ and Si NC with D = 3 nm [6j are also shown. Note that the 
experimental result for Si NCs [H] is higher than the bulk values. Our result 
for the NC with D = 3 nm is in good agreement with the experimental NC 
value. Our calculations lie even above the experimental ones. However, we 
think that a better agreement will be obtained when the experiments are held 
at a broad range of laser wavelengths. As a result, the calculations show that Si 
NCs access large values of the ratio n2//3A, particularly at small NC volumes. 
This feature should be taken into account in assessing Si NCs for nonlinear 
device applications especially for optical switching systems p[22] . 

In summary, we have investigated the wavelength and size dependence of the 
bound-charge third-order optical nonlinearities in Si NCs. It is observed that 
both n2 and /3 are enhanced with the decreasing of the NC size, giving values 
greater than their respective bulk values. Finally, optical switching parameter 
is assessed based on the numerical results. Si NCs enhance this parameter 
with respect to bulk Si. 
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